Dynamically generated dimension reduction and crossover in a spin orbital model 



o 
o 



Theja N. De Silva a ' b , Michael Ma a ' c , and Fu Chun Zhang a ' d 
a Department of Physics, University of Cincinnati, Cincinnati Ohio 45221 

b Department of physics, University of Ruhuna, Matara, Sri Lanka 
c Department of Physics, Hong Kong University of Science and Technology 
Department of physics, University of Hong Kong, Hong Kong 

We study a spin orbital model in which the spin-spin interaction couples linearly to the orbital 
isospin. Fluctuations drive the transition from paramagnetic state to C type ordered state into a 
strongly first order one, as observed in V2O3. At T = 0, there is a ferro-orbital-C-spin to ferro- 
orbital-G-spin transition. Close to the transition point, the system shows dynamically generated 
dimension reduction and crossover, resulting in one or more spin reentrant transitions. 
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Recently there has been growing interest in the effects 
of orbital degeneracy on the physics of transition metal 
oxides. The effective spin Hamiltonian of the insulating 
phase of such systems may depend crucially on orbital 
short-ranged and long-ranged correlations. As a result, 
magnetic ordering can become anomalous or may even be 
suppressed altogether. At the same time, orbital physics 
is also affected by spin fluctuations and correlations. The 
interplay between spin and orbital degrees of freedom is 
fundamental to much of the physics of transition metal 
oxides HQ. In this paper we investigate certain aspects 
of this interplay, emphasizing on the effects of thermal 
and quantum spin fluctuations on orbital ordering and 
the effects of orbital ordering on spin physics. Our most 
interesting result is that the system can exhibit dimen- 
sion reduction and dimension crossover of spin physics as 
a function of temperature due to orbital ordering. 

In the insulating phase of transition metal oxides, the 
dominating energy scales for the transition metal ions 
are the on-site Coulomb repulsion, Hund's rule cou- 
pling, and the crystal field due to the surrounding oxy- 
gen ions. Neglecting weak spin-orbit effect, the general 
spin Hamiltonian with two-fold degenerate orbitals (rep- 
resented by pseudospin r = 1/2) is of the form [jj , 
H = (JijSi ■ Sj + Kij), where Jy and are func- 
tions of Tj, and Tj. This Hamiltonian has global SU(2) 
invariance in spin space and a lower and discrete rota- 
tional symmetry in r space 0, HJ 0, 0, 0]- Out of the 
general class of such Hamiltonians, we will focus on those 
on a cubic lattice of the form 



H = J Si • Sj - K s * • s o(n 



na + Tj 



(1) 



Here the unit vectors fiy = fii, ri2, 83 for i,j nearest 
neighbor in the x, y, z directions respectively. In this 
model, the interplay between spins and orbitals arises 
from the second term which is linear in r. This linear 
term will be present provided the two eigenvalues of the 
hopping matrix are different, while the fij's will depend 
on how the two degenerate orbitals transform under lat- 
tice rotations. For specificity, we take the 's to be 
unit vectors in the x — z plane, with 113 = z, while Hi 



and S2 are rotated from S3 by 120° and 240°. More gen- 
erally, there are also quadratic in r terms, which we as- 
sume to be weak compared to the two terms kept. This 
implies that any orbital ordering in the system will be 
due to spin-orbital coupling rather than Jahn- Teller ef- 
fect. Assuming this is the case, and with the choice of 
rii 's above, this spin-orbital Hamiltonian can serve as 
a model for one electron or hole per site in the doubly 
degenerate e g levels of cubic peroskites Q as well as a 
possible model for V2O3 [Sj. For the latter, each site on 
the cubic lattice is the topological equivalence of a ver- 
tical pair of sites on the corundum lattice of V2O3. The 
coupling Jo depends strongly on and decreases with the 
Hund's coupling, while K is only weakly dependent on 
it. We consider Jo > 0, and with an appropriate defini- 
tion of r, we also have K > 0. The calculation shown 
here will be for S = 2, the value of S for the V2O3 bond 
model, but the results are qualitatively the same for other 
S. The results are also applicable to other lattices and 
other choices of the fij's. 

While Jo favors conventional (G-type) AF magnetic 
correlations so that nearest neighbor spins are all AF cor- 
related, K favors, along with orbital ordering, anomalous 
magnetic correlations break the cubic lattice rotational 
symmetry, for example C-type ordering with AF correla- 
tion in ab plane and FM correlation in c direction. In this 
paper, we investigate the phase diagram of this model in 
the temperature T and Jq/K plane. Our main results 
are: 1) At low T, orbital ordering gives rise to effective 
dimension reduction of the spin physics for Jq/K close to 
2. 2) The weakening of orbital ordering with increasing 
T can lead to dimension crossover from 2D to 3D and 
vice versa. 3) The dimension crossover effect together 
with thermal fluctuation effects on the spins can lead to 
an order by disorder mechanism and one or more reen- 
trant transitions. 4) The strongly first order nature of 
the magnetic transition in V2O3 is explained 10]. The 
underlying physics behind these results are orbital or- 
dering coupled with spin short-ranged correlations and 
quantum fluctuations. 

Within the context of the bond model for V2O3, our 
model has been studied by Joshi et. al. using a single-site 
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mean field theory. In order to include short range spin 
correlation and quantum fluctuations, we use a modified 
mean field approach to decouple the spin orbital Hamilto- 
nian. We begin with the Feynman-Hellman theorem : 

F < F + (H - H ) (2) 

Here F is the true free energy of the system and H is 
the actual Hamiltonian given in Eq. 1. Hq is a vari- 
ational Hamiltonian and () is the thermal expectation 
value with respect to Hq. Fo is the free energy of the 
system with Hamiltonian Hq. We take our variational 
Hamiltonian as, 

H Q =H S + H T (3) 

z x,y 
(ij) (ij) 

i 

where, Jj_ = Jn — 2Kt, and Jn = J + Kt, with t > 0. 
In H$ above, the first sum is for nearest neighbor pairs 
along the z direction and the second sum is for pairs on 
the same xy plane, t and A are variational parameters. 
Our choice of Hq is based on the expectation that orbital 
ordering will be ferro-orbital. 

Minimization of Feynman-Hellman free energy with re- 
spect to the variational parameters gives two self consis- 
tent equations, A = AB and t = (t z ) = ta,vii(2K/3A). 
AB = B± - B h with B ± = (Si • Sj) z and B\\ = 
(Si ■ Sj) x as out of plane and in-plane nearest neigh- 
bor spin spin correlation obtained from Hg. 

Note that while non zero value of t signifies long range 
orbital order, non zero value of B only signifies short 
range spin correlations. If B±_ ^ -By , spin correlations 
will be different from isotropic G type antiferromagnet. 
In particular, if Bu < 0, and B± > 0, that would corre- 
spond to C type magnetic correlations. Our decoupling 
scheme allows us to study the effects of short-ranged spin- 
spin correlations, but not short-ranged orbital correla- 
tions. However, we expect spin fluctuations to be dom- 
inant because of its continuous symmetry. Joint spin- 
orbital correlations are also ignored, but we expect these 
to be weak compared to spin fluctuations far away from 
the SU(A) limit \vi]. 

Note that Jii > 0, but J± can be either positive or 
negative in Hs, depending on the value of orbital or- 
der parameter t. Hs is an Heisenberg Hamiltonian with 
spatial anisotropy. In order to include short-ranged cor- 
relations and quantum fluctuations, we use renormal- 
ized spin wave theory (RSWT) [H Q Hljj . which un- 
like traditional spin wave theory (SWT) is applicable to 
both magnetically ordered and disordered phase. For the 
Heisenberg antiferromagnet on the square lattice, Hirsch 
and Tang [J^ have shown that this method can pro- 
vide quantitatively accurate results. In RSWT, magnon- 
magnon interactions are approximated by introducing a 



constraint that the total staggered magnetization be zero; 
M = J2ieA &i _ J2jeB Sj ~ 0, or equivalently that the 
average number of spin waves per site is S. This con- 
straint can be implemented by introducing into Hs a 
Lagrange multiplier A: 

z x,y 

H S = J ± J2 ' Sj + J\\ Si ■ S i ~ XM ( 4 ) 

(ij) <*i> 

The modified Hs (Eq. 4) is then solved using usual spin 
wave theory by expanding to quadratic order in Holstcin- 
Primakoff bosons. 

Classically, the spins will order as C-type and G-type 
for the J± < and J± > cases respectively. Thus, 
the sublattice designation will differ in the two cases and 
the spin wave calculation must be done separately. Let 
us define Q — ^ an d ^ and uj^ , magnon energies for 
Q > and Q < cases respectively. We see that Q is 
a dimensionless measure of the effective spin-spin cou- 
pling anisotropy. Note that — 2 < Q < 1. The chemical 
potential \i is obtained from the constraint equation, 

„ 1 f d d k ( 1 \\ , 

where (5 is the inverse temperature and u>^ are 
given by, ut = QJnS-\/ y? — jk 2 and uj^ — 

6J\\S^ (fi — ^7fe±) 2 — (|7fc||) 2 with effective magnon 
chemical potential, fi = |(2 + \Q\ + j^j-)- Ik = 
i(cosfcj. + cos k y + Qcosk z ), jk± — cosfc z and 7j.|| = 
^(cosfcj, + cosfc y ). F + (fi,k) = (i([i 2 — ^ 2 an( i 
F-fa,k) = ( / x+f 7fe± )(( M +| 7fc± ) 2 -(| 7fc1| ) 2 )-V2 should 
be used for Q > 0, < respectively. 

After solving for fj,, AB can be calculated, and the 
self consistent equation for t can be solved using itera- 
tive scheme. When there are more than one solutions 
for t, we compare their free energies to choose the sta- 
ble solution for each temperature. One anomaly of the 
RSWT approach is that B± and Bn need to be calcu- 
lated to one order of l/S higher than the free energy 
to ensure the correct sign for spin correlation at high 
temperatures 0, This difference is not significant 

for the self-consistent solutions close to and below the 
temperature of the first phase transition in our model. 

In FIG. 3, we show the resulting phase diagram for 
the case of S = 2, the value of S for the V2O3 bond 
model. In what follows, P stands for para (i.e. disor- 
dered), F stands for ferro ordering, C stands for C type 
ordering, G stands for G-type ordering, and O and S refer 
to orbital and spin respectively. Five phases are possible 
in the model. These phases are POPS , FOCS, FOPS 
(may be GS or CS short-ranged correlations), POGS 
(isotropic AFLRO), FOGS (anisotropic AFLRO). The 
phase consistent with the magnetic and orbital ordering 
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observed in V 2 3 is the FOCS phase HHE1011. The 

phase transitions the system undergoes as the tempera- 
ture is lowered depend on the parameter Jq/K and can 
be grouped into 6 regimes discussed below. In all cases, 
the orbital transition is first order while the spin transi- 
tions are second order unless accompanied by the orbital 
transition. 
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FIG. 1: (a) Phase diagram of spin-2 doubly degenerated spin 
orbital model, (b) Blow up of the region close to Jo/K = 2. 
Regimes III, IV, and V described in the text correspond to 
values of Jo/K between AB, BC, and CD respectively. Dotted 
line gives T2, the dimension reduction temperature. To its left 
(right), the spin coupling anisotropy parameter Q < 0(> 0). 

For regime I (0 < $ < 0.506), the regime relevant 
for V2O3, the transition is a direct one from the disor- 
dered phase into the FOCS phase with both spins and 
orbitals ordered. Experiments showed this transition to 
be strongly first order. While single site mean field the- 
ory (SSMFT) for our Hamiltonian is able to obtain a first 
order transition in this parameter regime, it was only 
weakly first order. By including short-ranged correla- 
tions in the present theory, the POPS phase is stabilized, 
and the first order nature of the transition is significantly 
enhanced. Entropy jump calculated from the free energy 
derivative, is larger than Iks per site throughout regime 
I, with e.g. a value of 1.28fcs for Jq/K = 0.25, which 
compares favorably with the experimental value of lAks 
(each site on the bond model corresponds to a pair of 
V ions). We should note that the mechanism for large 
entropy jump here is quite different from the usual fluc- 



tuation driven first order transition where the dominant 
fluctuations are from modes close to the mean field free 
energy minimum. For this system, however, the dom- 
inant fluctuations are G-type magnetic ones, and very 
far from the C-type magnetic ordering. Another issue in 
V2O3 is why spins and orbitals order at the same temper- 
ature. In principle, it is possible to have a spin Peierls 
transition driven by FO ordering instead of a FOCS tran- 
sition 6J. This issue cannot be addressed by SSMFT, 
but can be using our modified MFT, which shows that 
for Jq/K appropriate for V2O3, there is no FO driven 
spin Peierls phase. However, we will see below that such 
a phase can indeed occur in other parameter regimes of 
our model, but it will necessarily be preceded by a POGS 
ordering at higher temperature. Thus, a phenomenolog- 
ical explanation of the concurrence of orbital and CS or- 
dering in V2O3 is that there is no GS ordering at a higher 
temperature. 

Regimes II and VI (0.506 < # < 1.976 and 2.013 < £ 
respectively) show two phase transitions. The system 
first undergoes an isotropic (Q = 1) POGS ordering. 
Then at a lower temperature, the spin-orbital coupling 
causes a first order FO transition that converts the spin 
ordering to CS (Q < 0) in regime II and aniostropic GS 
(Q > 0, 7^ 1) in regimes VI. Rather more interesting, 
however, are regimes III, IV, and V, corresponding to 
1.976 < # < 1-987, 1.987 < ^ < 2, and 2 < ^ < 2.013 
respectively. These regimes show multiple transitions, 
including reentrance. These transitions are consequences 
of effective dimension reduction and dimension crossover 
in spin physics caused by orbital ordering. 

To see this, let us first consider T = 0, where within our 
MFT the orbital is always fully ordered. As a result, the 
effective spin Hamiltonian from our decoupling scheme 
parameterized by Q changes continuously to smaller pos- 
itive Q and then eventually to negative Q as Jo /K is de- 
creased. At Jq/K = 2, Q = and the spins on different 
planes become decoupled, i.e. the spin Hamiltonian is 
that of a 2D Heisenberg antiferromagnet. This is the or- 
bital driven dimension reduction effect at T — 0. Current 
wisdom is that the ground state of the 2D Heisenberg an- 
itferromagnetic Hamiltonian is ordered at T = in 2D 
even for S — 1/2. Thus there is spin LRO (anisotropic 
GS or CS) for all values of Jq/K at T = 0. At finite 
temperature, thermal fluctuations will weaken both the 
spin ordering and the orbital ordering, with the latter giv- 
ing rise to effective temperature-dependent spin Hamilto- 
nian. The spin physics is best understood by considering 
how Q(t) changes with T together with the dependence 
of the spin transition temperature T C {Q) on Q. For small 
\Q\,T C (Q) ~ \Q\ 9 ■ RSWT gives 6 = 1/2. Because or- 
bital has a discrete symmetry, its order parameter t, and 
hence Q changes exponentially slowly at low T, and the 
physics is dominated by thermal disordering of the spin. 
At higher T, the reduction in t becomes significant, and 
the corresponding change in Q can give rise to dimension 
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crossover in spin physics. These features are shown in 
Fig. 3b and discussed below, where Qo denotes the value 
of Q for t = 1, i.e. at T = 0. 

We first consider right at the decoupling point J /K = 
2 (point C in Fig. 3b), so that Q a = at T = 0. As 
T increases, t decreases and Q becomes increasingly pos- 
itive, the planes become increasingly coupled, implying 
a crossover from 2D to 3D. However, this crossover is 
exponentially slow at low T, and since Heisenbcrg spins 
cannot order at any T > in 2D, the spin LRO is im- 
mediately destroyed at infinitesimal T. As temperature 
increases, Q becomes large enough that T C (Q) exceeds T, 
and there is a reentrant transition into an anisotropic GS 
phase. The restoration of spin LRO due to temperature 
induced dimension crossover can be viewed as a novel 
kind of order by disorder mechanism. For Qo > but 
small (regime V), the physics is basically the same with 
one difference. Since now T c (Qo) > 0, the GS order is 
stable at low T but will disorder for T > T c (Qo). 

The behavior is even richer for Qo < but small 
(regimes III and IV). Now as T increases, Q gets first 
less negative, becomes at some temperature T2, then 
becomes positive. That is we have dimension crossovers 
first from anisotropic 3D to 2D and then back to 
anisotropic 3D as T increases. Correspondingly, the spins 
first undergoes a transition from CS LRO to CS short 
ranged order at T w T c (<3o)- The interplane ferrromag- 
netic correlation continues to decrease as T increases, 
crossing over into GS short-ranged order for T > T^. 
In regime IV, there is yet another reentrant transition 
into anisotropic GS LRO. Throughout regimes III - V, 
the PS phases has short-ranged spin correlations that 
are spatially anisotropic and so break the lattice rota- 
tional symmetry. In effect, these are orbital driven spin 
Peiriels phases. Eventually, in all these regimes, the sys- 
tem switches back in a first order jump back into isotropic 
GS ordering when the orbital becomes disordered. Point 
X is where this transition coincides with the dimension 
reduction temperature T2, so amazingly there is a jump 
directly from isotropic 3D behavior just above the tran- 
sition to exactly 2D just below. 

In summary, we have investigated the problem of the 
interplay between spins and orbitals in transition metal 
oxides concentrating on the competition between spin- 
spin interactions and spin-orbital coupling. In addition to 
illuminating the phase transition properties of magnetic 
ordering in insulating V2O3, our model shows a mech- 
anism for dynamically generated dimensional reduction 
and dimension crossover. Although the results presented 
are for S = 2, the same qualitative behavior will hold 
for other S. Also, while our calculations are restricted to 
the Hamiltonian (Eq. 1), these effects will be present in 
other spin-orbital models as long as orbital ordering re- 



sults in vanishing spin-spin coupling in one or more spa- 
tial directions. Since these dimensional reduction and 
crossover effects are present only close to the decoupling 
point, to observe them one would need to find systems 
with the appropriate Hund's coupling so as to produce 
the proper Jo/K range. For S = 2, there is the addi- 
tional problem that this range is very narrow. Larger 
range will occur for smaller S. Therefore it will be inter- 
esting to search for effective 5=1/2 peroskite transition 
metal oxides with double orbital degeneracy and weak 
Jahn- Teller coupling. 
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